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The possible infrared behaviour of the gluon propagator is stud- 
ied analytically, using the Schwinger- Dyson equations, in both the 
axial and the Landau gauge. The possibility of a gluon propagator 
less singular than 1/k 2 when k 2 — > is investigated and found to 
be inconsistent, despite claims to the contrary, whereas an infrared 
enhanced one is consistent. The implications for confinement are 
discussed. 
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1. Introduction 



The gluon propagator A^k) is gauge dependent and as such is not experimentally 
observable. However it's infrared behaviour has important implications for quark con- 
finement. It can be shown, that a gluon propagator, which is as singular as 1/k 4 when 
k 2 — > indicates that the interquark potential rises linearly with the separation. More 
formally, West [|TJ proved, that if, in any gauge, A^ v is as singular as l/k A then the Wilson 
operator satisfies an area law, often regarded as a signal for confinement. Hence, quarks 
are confined through gluon interaction. 

Another sufficient condition for confinement is, that a propagator of a coloured state 
should not have any singularities on the real, positive /c 2 -axis ||. So, if gluons are confined, 
then they cannot propagate on-shell and A^ must be less singular than 1/k 2 when k 2 — > 0. 
Such a behaviour of the gluon propagator was first assumed by Landshoff and Nachtmann 
|| on purely phenomenological grounds, being needed to reproduce experiment with their 
model of the pomeron. 

Clearly (Fig. 1) the gluon propagator cannot both be more singular and less singular 
than 1/k 2 as k 2 — > 0, but which is correct ? The Schwinger- Dyson equations provide the 
natural starting point for a non-perturbative investigation of this infrared behaviour of 
the gluon propagator. Extensive work has been previously performed in both the axial 
gauge ff]-[0] and the Landau gauge |8|-p0|. (For a comprehensive review see Roberts 
and Williams |pT|.) A solution as singular as l/k A has been shown to exist in both gauges 
H-H and KJ-Rnj, whereas a confined solution for the gluon propagator, i.e. less singular 
than 1/k 2 , has only been claimed to exist in the axial gauge J7|. The purpose of this 
paper is to explore why these two different behaviours have been found. Fortunately, 
in studying just the infrared behaviour, there is no need to solve the Schwinger-Dyson 
equation at all momenta. It is this that greatly simplifies our discussion and allows an 
analytic treatment. 

In Sect. 2 we briefly describe the Schwinger-Dyson equation for the gluon propagator. 
The axial gauge studies are reviewed in Sect. 3 and the possible, self-consistent solutions 
for the infrared behaviour of the gluon propagator are reproduced analytically. In Sect. 4 
we repeat the discussion for the Landau gauge and find, that a propagator less singular 
than 1/k 2 when k 2 — > is not a solution of the Schwinger-Dyson equation. In Sect. 5 we 
discuss the differing forms of the Schwinger-Dyson equations used to deduce these results. 
In Sect. 6 we state our conclusions. 
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Figure 1: Possible behaviour of the gluon propagator A(k 2 ), which is the coefficient of 
the g^ u or 5^ u component of /S.^ik). 

(a) confining gluon, A ~ (k 2 )~ 2 , 

(b) confined gluon, A ~ (k 2 )~ c with c very small, 

(c) infrared vanishing gluon A ~ k 2 . 

All are matched to the perturbative behaviour for k larger than a few GeV. 
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2. Schwinger-Dyson equation for the gluon propagator 

The Schwinger-Dyson equations are coupled integral equations which inter-relate the 
Green's functions of a field theory. Since they build an infinite tower of coupled equa- 
tions, approximations and truncations are necessary to solve them. The Schwinger-Dyson 
equation for the gluon propagator yields a relation for A M „ in terms of the full 3 and 
4-point vertex functions, and 1^ , the quark and the ghost propagators and cou- 
plings. The equation is displayed diagrammatically in Fig. 2. Here we only consider a 
pure gauge theory, i.e. a world without quarks. This is reasonable since we expect it is 
the non-Abelian nature of QCD which is responsible for confinement. 




Figure 2: The Schwinger-Dyson equation for the gluon propagator 

Here the broken line represents the ghost propagator. The symmetry factors 1/2 and 1/6 
and a negative sign for every ghost and fermion loop arise from the usual Feynman rules. 
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3. Axial Gauge Studies 



In the axial gauge the gluon propagator is transverse to the gauge vector n M , so 

Axial gauge formalism : n fl A afl = 

Studies of the axial gauge Schwinger- Dyson equation have the advantage that ghost fields 
are absent and the four-gluon vertex terms, Fig. 2, may be projected out of the Schwinger- 
Dyson equation. However they have the drawback that the gluon propagator depends not 
only on p 2 , but also on the unphysical gauge parameter 7, defined as 

(n ■ p) 2 
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n 2p2 



and in general must depend on two scalar functions, F and H 



1 



A„(p 2 , 7 ) = ~~ 2 [F(p 2 , 1 )M flu + H(p 2 , 1 )N^\ , (1) 



with the tensors given by: 



Mfu, = g^ u h n 



n ■ p (n • p) 2 



The free propagator is obtained by substituting F = 1 and H = 0. 

In all previous axial gauge studies it has been assumed that any infrared singular part 
of the propagator has the same tensor structure as the free one (though importantly this 
contradicts the results of West [13|]). Thus for p 2 — > it is assumed that 

A, u (p 2 , 1 ) = -^F(p 2 , 1 )M, u . (2) 

The gluon vacuum polarization tensor Yl\^(p 2 ,j) is defined by 

A„ = g\ - ^ . 
n ■ p 
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Projecting the integral equation with n^n^/n 2 the loops involving the four gluon vertex 
give an identically zero contribution because of the tensor structure of the bare 4-gluon 
vertex and the fact that the gluon propagator is transverse to the axial gauge vector, that 
is 

Uu^T = = A^n„ 
Thus the relevant part of the Schwinger- Dyson equation of Fig. 2 becomes: 

- 9 ^j-0^K^(p,k,-k,- P )A^(k) (3) 

where q = p — k, the last term is the tadpole contribution and all colour indices are 
implicitly included in the vertices. Once the full 3-gluon vertex is known, we have a 
closed equation for the gluon vacuum polarization IL^. 

The vertex is constrained by the Slavnov- Taylor identity in terms of this vacuum 
polarization : 

q„T pup (p, q, k) = W p (k) - W p (p) . (4) 

Separating T^ up into transverse and longitudinal part, where the transverse part is defined 
to vanish when contracted with any external momentum, the Slavnov- Taylor identity 
exactly determines the longitudinal part [T3] if it is to be free of kinematic singularities. 



Thus 



rL P (p,fc>?) = g^' !>! ' /: 



^\F(p\y) F(g2, 7 ), 



+ cyclic permutations . (5) 



This longitudinal part is responsible for the dominant ultraviolet structure of the 
vertex. Moreover, it is assumed, that it entirely embodies the infrared behaviour, and so 
the transverse part can be neglected. This assumption is motivated by the fact that the 
transverse part (as defined) vanishes, when the external momenta approach zero. 

Using the explicit expressions, Eq. (2) for A, Eqs. (4,5) for V and multiplying with 
n^riy/n 2 , we find in Euclidean space: 
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P 



F(p 2 



(1-7) 



2/1 \ I 9 2 C A f d*k n-(k-q) x A 
-V (1 - 7) + / j^y ~ 2 \o) (k) A (0) (g) 



fc ■ n 



2 J l^TTj" 1 

, F(q 2 ) - F(fe 2 ) 
J fc 2 - g 2 

- F(p 2 ) F{q 2 



(5 pa k ■ q - k p q a ) 
Pa(p + k) p 



-q ■ n 



p 2 - k 2 F{p 2 ) 
S pr7 F(k 2 ) - - ^ ' — ' v " ' [b pa k ■ q - k p q a ) 



n</ 2 ) ~ i'Xk 2 ) 
k 2 — q 2 

+ F( q 2 )-F(p 2 )F(k 2 ) 

+ p 2_ q 2 F{p2 f P K<l + P)v 



| g 2 C A I d 4 k F(k 2 ) ( 2 | Pn 2 



(2vr) 4 k 2 



(n ■ k)' 



(6) 



This is the equation first found by Baker, Ball and Zachariasen [|j] who studied its solution 
numerically. They came to the conclusion that the only consistent infrared behaviour for 
the function F(p 2 ) is 



F(p 2 ) oc — as p 
p 1 







and that this is independent of 7 as a numerical approximation. 

Schoenmaker |J simplified the BBZ equation further by exactly setting 7 = 0. Doing 
this the contribution of the tadpole diagram vanishes. Moreover, approximating F(q 2 ) by 
F(p 2 + k 2 ), which should be exact in the infrared limit, allows the angular integrals to be 
performed analytically. Consequently, Schoenmaker finds the following simpler equation: 



P 



F(p 2 ) 
'a 



1 



fa 

32tt 5 



dk 2 



k A 5 k 2 2 k 2 s 
12/7 + 2f ~ 3p 2 - k 2 



Ft 



6 p 2 
3 k 2 



k 2 



3 p 2 — k 2 



POO 

+ / dk 2 

Jp 2 



( k 6 _lk A 
f \2Af ~ 4p* 
3 p 2 2 p 2 



lk 2 



Ap 2 



I — I 5 lp2 ) F I ( 7p2 I 2 P 
I Ap 2 12 8 k 2 2 UP 3p 2 -k 2 



A Ak 2 3p 2 -k 2 
2 \ 



F L 



(7) 



where 



F 2 
F 3 



F(p 2 + k 2 ) 

F(p 2 + k 2 ) - F(k 2 ) 
F(p 2 + k 2 )F{k 2 ) 
F(p 2 ) 



In general, this equation has a quadratic ultraviolet divergence, which would give a mass 
to the gluon. Such terms have to be subtracted to ensure the masslessness condition 
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lim IV = , i.e. =0 for p 2 ^0 , (8) 

is satisfied. This property can be derived generally from the Slavnov- Taylor identity and 
always has to hold. 

The complicated structure of the integral equation, Eq. (7), does not allow an exact 
analytic solution for the gluon renormalization function F(p 2 ) to be found and most 
previous studies (||, |K|, (4| and 0) solve the equation numerically. However, the 
possible asymptotic behaviour of F(p 2 ) for both small and large p 2 can be investigated 
analytically. 

We determine which infrared behaviour of F(p 2 ) can give a self-consistent solution to 
the integral equation by taking a trial input function F in (p 2 ) and substituting it into the 
right hand side of the equation. After performing the fc 2 -integration, we obtain an output 
function l/F out (p 2 ) to be compared to the reciprocal of the input function. To do this, 
the gluon renormalization function is approximated in the infrared region by a Laurent 
expansion in powers of p 2 and at large momenta by its bare form, i.e. 

En= an(p 2 /H 2 ) n+V for p 2 <^ 2 

nin = { , (9) 

1 for p 2 > ft 2 



where 



n=0 



to ensure continuity at p 2 = /i 2 . /i is the mass scale above which we assume perturbation 
theory applies, rj can be negative to allow for an infrared enhancement. Eq. (9) is a suffi- 
ciently general representation for finding the dominant self-consistent infrared behaviour. 
Of course, the true renormalization function is modulated by powers of logarithms of 
momentum, characteristic of a gauge theory. However, these do not qualitatively affect 
the dominant infrared behaviour and can be neglected. Indeed to make the presentation 
straightforward, we only need approximate F(p 2 ) by its dominant infrared power (p 2 ) v 
for p 2 < fi 2 to test whether consistency is possible and this is what we describe below. 
However, as we shall see, if t] is negative then potential mass terms arise and these have 
to be subtracted. Only in this case do higher terms in Eq. (9) play a role too and it is 
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necessary to consider other than the leading term in the low momentum input. Other- 
wise higher powers make no qualitative difference as we have checked. Consequently we 
present only the results with the lowest powers in the representation, Eq. (9). 

To illustrate the idea, let us take the trial infrared behaviour to be just 



' 2 \ V 
pi \ 



F(p 2 ) oc (i.e. a n = for n > 1) . (10) 

Note that the masslessness condition, Eq. (8), restricts r] to be less than 1. Furthermore we 
demand that in the high momentum region the solution of the integral equation matches 
the perturbative result, i.e. for p 2 — > oo, we have F(p 2 ) = 1, modulo logarithms. 

Taking 77 = —1, for example, i.e. 

1 2 
P 2 

in Schoenmaker's approximation, Eq. (7), gives 



FUP 2 ) = 



This violates the masslessness condition of Eq. (8) and so has to be mass renormalized. 
As explained above, now terms in F(p 2 ) of higher order in p 2 will generate a contribution 
to the right hand side of the equation making it possible to find a self-consistent solution 
by these cancelling the explicit factor of 1. Consequently, we can approximate Eq. (9) by 

A (/i 2 /V) + (p 2 //J> 2 ) if p 2 < H 2 
F in (p 2 ) = { , (11) 



if p 2 > p 2 



We then find, after mass renormalization: 



1 1 + 9 *c A 



96 /i 2 12 8/j 2 n [p 2 ) + 12 n I n 2 



(12) 



F out (p 2 ) 32vr2 

where A is the ultraviolet cut-off introduced to make the integrals finite. The ultraviolet 
divergent constant can be arranged to cancel the 1 and we find self-consistency modulo 
logarithms. It is this result that Schoenmaker found P supporting the earlier result of 
BBZ [f|. However, importantly, self-consistency requires A, Eq. (11), to be negative as 
also found by Schoenmaker. 
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More recently, Cudell and Ross [|7| have taken Schoenmaker's equation, Eq. (7), and 
investigated whether one can find self-consistency for a gluon renormalization function 
which is less singular than 1/fc 2 for k 2 —> 0, i.e. which corresponds to confined gluons. 
The trial input function they use in their investigation is 



F m {p 2 ) « ( P 2 y- c , 

where c is small and positive to ensure a massless gluon, Eq. (8) . Once more we want the 
integral equation for IT^ to agree with perturbation theory in the ultraviolet region, but 
Fin{p 2 ) oc (p 2 ) l ~ c grows for large momenta and hence spoils the ultraviolet behaviour. So 
to check whether this input function gives self-consistency in the infrared, we input the 
trial form : 

!(p 2 //i 2 ) 1_c if p 2 < /i 2 
(13) 
1 if p 2 > fi 2 

Inserting this into Eq. (7), we find, after mass renormalization : 



i 1 g 2 C; 
1 + 



Dl+ D Aff\ D jff\ D m\ 



,(14) 



F out (p 2 ) 327T 2 

where Fx,Fz and F% have been expanded for small p 2 and only the first few terms have 
been collected in this equation so that 
5 1 11 A 2 

D 1 = ; r 1 In — , 

12(1 -c) 3 12 /i 2 
1 1, /A 2 \ 



D 



2407 353c 3 1 l + 2c 7- 8c 7 + c 7 



3 



D 4 



1440 360 8c 24(5 - c) 12(4 - c) 12(3 - c) 12(2 - c) 12(1 - c) 

+ w^k + w^t) + ¥ 2 - c >*« - > - c )*p - c > + ¥ 2 - - 2c > ■ 

7 -9c 
24^ 



where ^ is the logarithmic derivative of the Gamma function. Again the 1 can be arranged 
to cancel with the constant term and the dominant infrared behaviour is indeed 



' o\ 1 — c 

/i 2 \ 



for p 2 -> (15) 



Fout(p 2 ) \P 2 J 

Hence a gluon propagator less singular than 1/p 2 for p 2 — > can be derived from Schoen- 
maker's equation as Cudell and Ross have found. Note once again that terms of higher 
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order in Eq. (9) do not qualitatively alter the result. Thus we see in the axial gauges 
that apparently both confined and confining solutions are possible for the gluon propaga- 
tor. However, the singular confining behaviour must be an artefact of the approximation 



that one of the gluon functions, H, vanishes, since West |Td| has shown that in a gauge 
with only positive norm-states, a singular gluon renormalization function is not possible. 
Moreover, the approximation of setting 7 = in the BBZ-equation, Eqs. (6,7), has been 



seriously questioned in Ref. |14|. It is therefore sensible to ask what is the behaviour in 



covariant gauges, to which we now turn. 
4. Landau Gauge Studies 

The advantage of Landau gauge studies is the much simpler structure of the gluon 
propagator, which is defined by : 

. G(p 2 ) ( p^p u \ 

A /- = ~ % —r- \9v --pir) ■ ( 16 ) 

However other problems arise and the following approximations have to be made : 



• In any covariant gauge, ghosts are necessary to keep the vacuum polarization trans- 
verse and hence are present in the Schwinger-Dyson equation of the gluon propaga- 



tor, Fig. 2. However, in all previous studies ||, [KJ the ghost loop diagram is only 
included in as much as it ensures the transversality of the gluon propagator, assum- 
ing that otherwise it does not affect the infrared behaviour of the propagator. This 
assumption is supported by the fact, that in a one- loop perturbative calculation the 
ghost loop makes a numerically small contribution to G(p 2 ). 

• The 4-gluon terms cannot be eliminated as in the axial gauge and are simply ne- 
glected. This can be regarded as a first step in a truncation of the Schwinger-Dyson 
equations. 



With these assumptions, we again find a closed integral equation for the gluon vacuum 
polarization, H^, once the full 3-gluon vertex is known. In the Landau gauge, the Slavnov- 
Taylor identity for the 3-gluon vertex involves the ghost self-energy, which is simply set to 
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zero, and the proper ghost-gluon vertex function G^ v . However, in the limit of vanishing 
ghost momentum the ghost-gluon vertex is approximately equal to the gluon propagator. 
With this simplification the Slavnov- Taylor identity has the same form as in the axial 
gauge and is given in Eq. (5). Once again neglecting the transverse part of the vertex, we 
obtain a closed integral equation : 



n, 



^ - ng - y / t Z(-p> k > ?) Aa "( k ) A %) -q) (17) 

where once again the colour indices are implicit and q = p — k. 
A scalar equation is obtained by projecting with 

P^ = -L(4p^-p 2 gn 

This projector has the advantage that the term in Eq. (16) that is quadratically 
divergent in 4-dimensions, does not contribute. Thus we find 

1 . 9 2 C A 1 



G{p 2 



1 + 



<fk 



96tt 4 p 2 
G(k 2 ) - G(p 2 ) G(q 2 
k 2 — p 2 G(p 2 



G(q 2 )A(k 2 y) + G(k y 9 [ q2) B(k 2 ,p 2 ) 



G{p 2 ) 

nn2 2\ i G(q 2 ) — G(k 2 ) „, ;2 2\ 

C{k ,p ) + — D(k ,p ) 

qz _ ^4 



where 
A(k 2 ,p 2 ) 

B(k 2 ,p 2 ) 
C(k 2 , P 2 ) 
D(k 2 ,p 2 ) 



Ig (fc-P) 2 _ 61 ( fc -P) + 16 ( fc "P) 3 _ 12 J_ + 22 P 2 
k 2 p 2 q 2 k 2 q 2 k 2 p 2 q 4 k 2 k 2 q 2 



42 



(k -pf 
k 2 q A 
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V 



k 2 q 4 



+36 



p 2 {k ■ p) 



T 



k 2 q 



k 2 q 4 



^k-pf 6 p 2 (k-p) 6 (fc-p) g p 



Fg 4 

,2 



fc 2 g 2 



£; 2 g 2 



12^! - ls ^ -p ) 2 + ls ( fc -p) 3 + 21 ( fc -p) _ 5 e! _ 10 ( fc -p) 2 + 9 f( k ■ P) 



q z 



p 2 q 2 



k 2 p 2 q 2 



q 



k 2 q 2 



k 2 q 2 



Brown and Pennington [10 studied this equation numerically and found 



G(p 2 ) =A^r for p 

p2 







to be a consistent solution. This result is in agreement with Mandelstam's study of the 
gluon propagator M. 
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Again approximating G(q 2 ) by G(p 2 + k 2 ) allows us to perform the angular integrals 
analytically in Eq. (18), giving: 
1 



1 + 



GV) 

9 2 C A 1 



487T 2 p 2 



dk 2 





G 2 



-1 



A; 2 A; 4 
10- + 6- + 



k 2 



75 
4 



39 A; 2 A; 4 
+ 4^ 



4 p 2 



P 



.p_ 
'k 2 



21k 2 k 4 k 6 \ „ 
-— — + 7— -3— ] +G 3 



27 



4 p 2 



pa 



+ 



+G 2 



dk A 



6 + 



> 2 - k 2 \ 

3p 2 ' 



P 



S P 

_ /29 
I 2 I T ' 4 A; 2 



11 A? 

4 p 2 



15 

IF 



where 



3 

2 + 4 A; 2 



Gi 
G 2 

G 3 



G, 



67p 2 _ 3p 4 * 
~8~P ~~ 8F 



G(p 2 + k 2 ) 

G(p 2 + k 2 ) - G(A; 2 ) 
G{k 2 )G{p 2 + A: 2 ) 
G(p 2 ) 



(19) 



Note that the integral equation has the usual ultraviolet divergences, but infrared 
divergences are also possible. The ultraviolet divergences can be handled in the standard 
way to give a renormalized function G#(p 2 ) — this will not be discussed here. However 
we have to make the potentially infrared divergent integrals finite in order to calculate the 
integrals^. The infrared regularization procedure proposed by Brown and Pennington 



10 is to use the plus prescription of the theory of distributions, which is defined as 



follows: 



ft. 

k 2 



k 2 



for oo > k > 



and in the neighbourhood of A; 2 = it is a distribution that satisfies: 

~2 Jl\ __ 



~S(k 2 ,p 2 )-S(0,p 2 ) 



p 2 u 2 
K 



fi 



2 k 



Simply taking 



G in (k 2 ) = A 



ft. 

k 2 



(20) 



1 These divergences do not arise in an axial gauge when 7 is set equal to zero as Schoenmaker does, 
Eq. (7). 
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as an input function once again leads to a mass term and higher terms in the expansion 
Eq. (8) are necessary. Then we do have the chance of finding self-consistency for a gluon 
propagator as singular as 1/fc 4 and hence confining quarks. With 

A <y /V )+ + (P 2 /V) if P 2 < l* 2 



G in (p 2 



if p 2 > /i 2 



(21) 



we find, after mass renormalization: 



G out (p 2 



9 2 C A 
48tt 2 



479 p 2 13 p 



fx 2 8/i 



81 _ 25 /A 2 ' 



4 



(22) 



The ultraviolet divergent constant can be arranged to cancel the 1 and, again, we find 



self-consistency. This is the result found numerically by Brown and Pennington |10j with 
a positive infrared enhancement to the gluon renormalization function, i.e. A > 0. 



Now we check, whether it is possible in the Landau gauge, to find the behaviour Cudell 
and Ross discovered using Schoenmaker's approximation in the axial gauge. With 

(p 2 //i 2 ) l ~ c ifp 2 </i 2 

1 if p 2 > /i 2 



G in (p 2 ) 



(23) 



we find, again after mass renormalization: 



G out (p 2 ) 
where 

D 2 
D 3 



Da 



1 + 



9 2 C A 
48tt 2 



'3 + 5 + 6c + 25 ln M 2N 



c 4 
3 



4(2 - 2c) 4 



In 



A 2 



13c 59 
+ 



32c 155 - 64c 
+ 



1971 29c 37 6 

IttT + ~Y + 20c + 2(1 - c) ' 4(2 - c) ' 8(3 - c) 
127 -49c 23 -11c 125 + 61c 55 + 6c 3 

+ ^7Z T - ^T, r-T - ^ =T + 



+ 



8(4 - c) 4(5 - c) 8(1 - 2c) 

-8(2 - c)*(-2c) - 8(2 - c)^(l) 

61 + 6c 
8(1 -2c) 



5(2 - 2c) 4(3 - 2c) 



(24) 
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Thus the dominant infrared behaviour is: 

— ! 

and self-consistency is spoiled by a negative sign, since c is small and positive. 
5. Confined Gluons 

A gluon propagator, which is less singular than 1/k 2 for k 2 — > 0, and hence describes 
confined gluons appears to be a self-consistent solution only of the axial gauge Schwinger- 
Dyson equation using Schoenmaker's approximate integral Eq. (7). In the Landau gauge 
this behaviour of the gluon propagator is not possible : a minus-sign spoils self-consistency. 
We should therefore comment on the origin of this crucial minus sign. 

Starting from BBZ's integral Eq. (6), there is no difference in sign between the two 
gauges. Eq. (6) is Bose-symmetric (as it should be) and can therefore be rewritten as : 

,r (l- 7 )=p 2 (l- 7 ) 



F{V 2 ) 

-^T / ^^^^)A^ (25) 



where 



K p „(k,p) = 5 pu F{q 2 ) - F{ ^ 2 _ F jf ] \ 5 pa k ■ q - k p q a ) 

F(k 2 ) - F{p 2 ) F{q 2 ) 
+ p 2 _ k 2 F(p2) IP + 

whereas, taking the starting equation of Schoenmaker's paper (Eq 3.5 of Ref. 0) we find: 
''" (l- 7 ) =^(1-7) 



F(p 2 ) 

~-/C A J ^ U ' ^ g)[~2fc • n K^Kp)} (26) 

where 

Z p(7 (k,q) = ( i ) 2 A^(k)A$ ) (q) . 

Schoenmaker formulates his equation in Minkowski space. Performing a Wick rotation to 
transform to Euclidean space by : d A ku —> id^kg, we find that Schoenmaker's equation, 
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Eq. (26), becomes : 



,2 



< r (l- 7 ) =^(1-7) 



F(p 2 ) 

^ I ^0^^^ (27) 

which differs from Eq. (25) by a crucial minus sign. 

We therefore see that in the axial gauge using BBZ's integral equation for the gluon 
propagator, and simplifying the angular dependence in the way Schoenmaker does in order 
to make an analytical discussion of the infrared behaviour of the propagator possible, 
yields an integral equation very similar to the one found by Brown and Pennington [|K^ 
in the Landau gauge. These equations lead to the correct perturbative behaviour at large 
momenta. In contrast a self-consistent solution of the gluon propagator less singular than 
1/k 2 for k 2 — > cannot be found in either gauge. Schoenmaker's own equation, which 
is the starting point for the study of Cudell and Ross |7| for instance, has an incorrect 
additional minus sign. This should have been heralded by the self-consistent enhanced 
gluon of Eq. (11) having a negative sign, using Schoenmaker's equation. In an axial gauge 
this sign should have been a little worrying for a wavefunction renormalization of a state 
with positive definite norm. 



6. Summary and Conclusion 

We have studied the Schwinger- Dyson equation of the gluon propagator to determine 
analytically the possible infrared solutions for the gluon renormalization function G(p 2 ). 
In both the axial and Landau gauges, one can find a self-consistent solution, which behaves 
as 1/k 2 for k 2 — > and hence a propagator which is as singular as 1/k 4 for k 2 — > 0. This 
form of the gluon propagator is consistent with area law behaviour of the Wilson loop, 
which is regarded as a signal for confinement 0. Numerical studies [|J, |10] have shown 



that a gluon propagator with such an enhanced behaviour in the infrared region that 
connects to the perturbative regime at a finite momentum (as indicated by experiment) 
can indeed be found as a self-consistent solution to the boson Schwinger-Dyson equation. 
Such a behaviour of the boson propagator has been shown to give quark propagators 



2 Though as remarked at the end of Sect. 4, axial gauge studies are seriously marred by the simplifying 
assumption that H = in Eq. (1). 
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with no physical poles [IS|. Furthermore, extending these non-perturbative methods to 
hadron physics, it has been found that a regularized, infrared singular gluon propagator 
together with the Schwinger- Dyson equation for the quark self-energy, gives rise to a good 
description of dynamical chiral symmetry breaking. For instance, one obtains values for 
quantities such as the pion decay constant that agree with experimental results [ |TT|1 . 

A gluon propagator which is less singular than 1/k 2 for k 2 — » 0, and hence describes 
confined gluons, cannot be found in either the axial or the Landau gauge. Solutions of this 
type have only been found using approximations to the gluon Schwinger-Dyson equation 
with an incorrect sign. Possible consequences for models of the pomeron are discussed 



elsewhere 19 



We should also mention the related work of the group of Stingl et al. ||15|| . They too 



start from an approximate, but larger, set of Schwinger-Dyson equations, which is then to 
be solved self-consistently. However, the method employed is completely different. The 
philosophy |T3[ is to obtain the solution of these equations as power series in the coupling, 



as in perturbation theory, and to include non-perturbative effects by letting each Green's 
function depend upon a spontaneously generated mass scale, b(g 2 ). The gluon propagator 
is assumed to be of the form (see Fig. 1) 

representing confined gluons. This grossly violates the masslessness condition of Eq. (8). 
In general, gluon masses can only arise in 4-dimensions if the vertex functions have singu- 
larities themselves corresponding to coloured massless scalar states, otherwise the Slavnov- 
Taylor identities sufficiently constrain the vertex functions to require the inverse of the 
gluon propagator to vanish at k — > []. Not only do the vertices of Stingl et al. have 
these massless singularities but self-consistency can only be found if the 3-gluon vertex is 
complex, when conventional understanding of its singularity structure would lead us to 
expect it to be real for momenta which in Minkowski space are spacelike. Subsequently 
Hawes, Roberts and Williams Jl6| and Alkofer and Bender [I7J] have shown that with this 
gluon propagator, Eq. (28), the solution of the quark Schwinger-Dyson equation does not 
describe a confined particle. They therefore also conclude that the full gluon propagator 
in QCD cannot vanish in the infrared region. 

3 This remark equally applies to the work of Zwanziger |20) and lattice studies by pl| . 
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To summarise: 

At first sight there appears to be a distinction between a confining and a confined gluon. 
A confining gluon is one whose interactions lead to quark confinement. A(k 2 ) ~ l/k 4 
behaviour is of this confining type. In contrast, it is sometimes argued that A(k 2 ) must 
be less singular than l/k 2 to ensure that gluons themselves do not propagate over large 
distances. However, whether gluons are confining or confined are not real alternatives. 
Gluons must be both. They confine quarks by having very strong long range interactions. 
They themselves are confined by not having a Lehman representation that any physical 
asymptotic state must have fll] • 

While infrared singular gluons satisfy both criteria, softened gluons though confined, 
do not generate quark confinement or dynamical chiral symmetry breaking, which are 
features of our world. Remarkably, a study of the field equations of QCD reveals this 
theory naturally exhibits these aspects with an infrared enhanced gluon propagator. 
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